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ABSTRACT 

A theoretical and experimental investigation to determine the 
lift coefficient of a flat plate planing surface covering a broad 
range of trim angles and aspect ratios is presented. 

The analysis is divided into two principal parts: first con- 
sidered is the modification of the two-dimensional potential flow 
solution for the effect of a finite span neglecting gravity; then, 
second, the effects of gravity on the finite span planing surface 
are discussed with the importance of Froude Number emphasized. 

The analysis made on the effect of a finite span neglecting 
gravity is accompanied by an experimental investigation in which 
angles of trim up to 30 degrees and wetted length-beam ratios up to 
about seven are explored. 

The analytical approach leads to the prediction of reasonable 
values for the lift coefficient, but to get better agreement with 
experimental data, an empirical correction factor is introduced. 

The investigation, conducted at the Experimental Towing Tank, 
Stevens Institute of Technology, Hoboken, N. J., is the thesis sub- 
mitted by the author in partial fulfillment of the requirements for 


the degree of Master of Science. 
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Geometric angle of trim 

Wetted length-beam ratio = 

Speed coefficient (Froude Number) V/vgb 
Acceleration of gravity 

Beam of planing surface or span of airfoil 


Velocity of the free stream 


— 25Ы 
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Lift coefficient, 5 
5 PV DL 


Lift force 

Fluid density 

Wetted length of planing surface or chord of airfoil 
e 

Aspect ratio, x Or = 

Projected wing area, | dy dx 

Normal force S 

Ambient flow angle of attack 

Induced angle of attack 

Spray thickness 

Induced downwash velocity infinitely far aft 

Virtual mass rate of flow, VeVol. 

Kinetic energy 

Change in kinetic energy 


Wetted length measured from trailing edge to stagnation 
line on planing surface 
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Empirical wave rise factor 


Wetted length-beam ratio neglecting wave rise 


Empirical correction factor 


Notation Used in the Section on the Effect of Gravity 
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Velocity potential 

Free stream velocity 

Coordinates of point (x,y) on a planing surface 
Coordinates of a point in space 


Wave number E, 
U 
Coefficient introduced by Maruo to obtain desired value 
of certain integrals 
Pressure at any point on planing surface 
Wave direction angle 
General wave number 
Downwash as used by Maruo 


A function of Froude Number 


Buoyancy force (Archimedian Force) 


Notation in Appendix I 


Geometric angle of attack 
Ambient flow angle of attack 
Induced angle of attack 
Gomplex potential 

Velocity potential 

Stream function 


Complex variable, z = x + iy 





T Kinetic energy 
Š Surface 


V Transverse flow velocity 


Notation in Appendix II 





R Force 

R Drag force 

X 

Q Volume flow in spray jet, Q = Vg 
ol Angle of trim 

Rç Spray resistance 

D Total lift force, dynamic plus Archimedian 
R Frictional Force 

tp 

R. Wave resistance 
y A function of Froude Number, = 1 5 

2 2 2 EE 

ET Froude Number squared, = 
M Moment 

С, Lift coefficient (C,) 
La Center of pressure location from trailing edge 


P S, R Vector forces 


h Depth of fluid from free surface 





THE LIFT COEFFICIENT OF FLAT PLAN ING SURFACES 





INTRODUCTION 


Hydrodynamic phenomena arising from planing on the surface of 
water have aroused considerable interest in seaplane hull, float, and 
hydroski enthusiasts. In the past the need for immediate design data 
for floats or hydroskis has led, based on many specific experimental 
investigations, to the development of numerous empirical formulae to 
describe the flat plate lift coefficient as it varies with other рага- 
meters such as angle of trim, Y, wetted length-beam ratio, A, and 
Froude Number, V/Vgb. While these empirical formulae agree very well 
with experimental results within the range of vaiues for which they 
were developed, the validity of such formulae outside the range of data 
for which they were developed becomes unreliable or at least question-- 
able. 

The objective of this work is to develop on a rational basis, 
insofar as possible, a formula for the lift coefficient of a flat 
plate which planes on the surfase of water. At present, a rationally 
derived formula developed by H. Wagner, Ref. 1, covers the infinite 
span case, neglecting gravity, while the more recent work of L. I. Sedov, 
Ref. 2, and H, Maruo, Ref. 3, covers the infinite span including the 


effect of fluid gravity which involves wavemaking as well as  hydrostatic 





effects. Methods of correction for the effect of a finite span neg- 
lecting gravity have been developed by such writers as B. V. Korvin- 
Kroukovsky, Byrne Perry, W. Bollay, W. G. A, Perring and L. Johnston, 
W. Sottorf, L. Landweber, F. W. S. Locke, and С. L. Shuford, Jr., 
Ref. l-12, but the entire process, including the effects of trim, 
aspect ratio, and Froude Number, has not been integrated except in 

a preliminary form. 

This work involves an investigation of the methods and applic- 
ability of derivations previously made by the investigators of flat- 
plate planing. The problem of deriving a formula on a rational basis 
is treated in two distinct parts. First, the aspect ratio, trim angle 
effect, neglecting gravity, is considered; then, second, the effects of 
gravity including wavemaking are discussed. 

To extend the range of experimental data for comparison with re- 
sults developed analytically, an experimental investigation was made 
which is described herein. 

The investigation has been conducted by Lt. D. D. Farshing, Jr., 
U.SoN,, a Naval Postgraduate School student at the Experimental Towing 
Tank, Stevens Institute of Technology, Hoboken, New Jersey.  Acknow- 
ledgement is made to Professor B. V. Korvin-Kroukovsky of the Exper- 
imental Towing Tank for his helpful guidance and interest in this work. 
Gratitude is also expressed to the members of the Towing Tank staff who 
have willingly given their assistance and advice throughout the course 


of this investigation. 
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MODIFICATION OF THE TWO-DIMENSIONAL PLANING SOLUTION 


FOR THE EFFECT OF A FINITE SPAN IN A WEIGHTLESS FLUID 


Historical Background 


There have been several attempts to find closed form solutions 
to the finite span planing surface problen, but difficulty arises in 
that formulation of the problem to satisfy the boundary conditions 
introduces a non-linear integral equation. Solution of the integral 
equation is feasible only after many simplifying assumptions have been 
made. 

To circumvent the difficulty arising from the integral equation, 
another approach is offered. An ideal weightless fluid is assumed to 
exist. Based on the method of potential flow, a two-dimensional solu- 
tion is found, then modified for the effect of a finite span by a tech- 
nique similar to that used for correcting the two-dimensional airfoil 
to finite span. 

A brief but comprehensive review of the significant work in plan- 
ing up to 1919 has been presented by William Siler, Ref. 13. Charles 
L. Shuford, Jr., Ref. 12, has given a recapitulation of the results of 
many investigators who sought to find a simple formula for the three— 
dimensional lift coefficient of the planing flat plate. After his re- 
view of the attempts of other investigators to find simple formulas, 
Shuford proposes another method for predicting the lift coefficient. 
Shuford observes that experimental data indicates a non-linear rela- 


tionship between the angle of trim and the planing lift coefficient 





for the planing surface of finite span. To account for the non- 
linearity, he approaches the problem with an analogy to the low aspect 
ratio airfoil. Appendix II shows further justification for this ap- 
proach. To arrive at a final total lift formula, Shuford combines 
three components; the first, a linear term given as one-half the value 
of the flat plate airfoil lift coefficient based on the lifting line 
theory; the second, a non-linear, cross flow term which is shown to be 
proportional to siny; and the third, a suction component of lift 
existing at the leading edge of the airfoil but not present in the 
planing surface flow. 

Shuford's total lift formula becomes 


2 


T 
U RT (1 - Sind) + sin T cos Т 


T: 1+AR 
values of which agree well with experimental data. The analysis by 
Shuford, while т good agreement with the available experimental data, 
is not completely theoretical since the planing problem was not sclved, 
but, instead, the airfoil problem manipulated. 

In many of the papers concerned with modification of the two- 
dimensional solution for the effect of finite span, methods similar to 
those used for correcting the infinite span airfoil have been used with 
varying degrees of success. The technique of adopting similar methods 
seems to be a logical approach since the flow is assumed ideal (without 
friction) and weightless. 

Another successful approach is one made by Byrne Perry, Ref. 6. 
Perry has applied the low aspect ratio airfoil analysis of Е. Weinig, 


Ref. ll, to the low aspect ratio planing surface. Weinig solved the 





airfoil problem without recourse to the two-dimensional solution. He 
applied the momentum theorem to a control surface around the airfoil 
with virtual mass and induced angle of attack as independent variables. 
The virtual mass used by Weinig is that associated with the lifting 
line plus a wedge of fluid found by the projection of wetted length 
normal to the flow: 
Жа ve (See Fig. 1) 

Appendix I presents more in detail the concept of virtual mass and the 
work of Weinig insofar as it is understood. In order to justify the 
choice of virtual mass made by Weinig, the applicable work of M. M. 
Munk, Ref. 15, is also discussed. The formulation of the low aspect 
ratio airfoil problem and its solution by William Bollay, Ref. 16, is 
considered appropriate for discussion in Appendix I, and it is com- 
pared with the solution obtained by Weinig. 

The principal objection to the planing surface analysis by anal- 
ogy with the low aspect ratio airfoil of Weinig is that the method for 


obtaining the induced angle of attack is not clear. 


Ihe Two-Dimensional Solution 


H. Wagner, whose applicable work on planing is found in Ref. 1, 
obtained a two-dimensional solution to the planing problem. As pre- 
viously stated, an inviscid, weightless fluid is assumed to exist. 
Coordinate axes are fixed in the planing surface which has infinite 
span and infinite forward length. Fluid is allowed to flow toward 


the plate parallel to the x axis shown in Fig. 2. The mass of fluid 





of thickness $ contained between the free surface and the stagnation 
streamline as shown on Fig. 2 is deflected forward along the plate. 
Ihe fluid below the stagnation streamline continues aft and finally 
becomes part of the wake. With the aid of the mathematics of complex 
variables using conformal transformations, Wagner finds the normal 


force per unit width of the infinite span to be: 


Care ysin Y 
s ¿AY 1~cosv 
l + cosg- (1 - cost)£n er UU Bin 


The projection of this hydrodynamic force perpendicular to the direc- 


tion of flow gives the lift force per unit span. Then: 


hoe Loy? 27 sin?’ cosv 
"2 l-cos i 
1 + соѕт- (1 – соѕт) п 522931 тіп у 
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For the infinite span case, Y -w, and ©. = 0. 


Application of the Momentum Theorem - Method I 


Let us now take the infinite span solution of Wagner and modify 
it for the effect of finite span. 

To begin we must put a control surface around the finite span 
planing surface so that the momentum theorem can be applied. For the 
control surface let us choose a solid shape of infinite spanwise 
width. One face is infinitely far ahead of the planing surface; one 
infinitely far aft; one at infinite depth below; and the remaining 


face at the free surface, a,b,c,d,e,f of Fig. 2. 





There is no momentum flux through the side and lower faces of 
the control surface, so we can focus our attention of the fore and aft 
faces and that part of the free surface face which cuts through the 
forward spray. 

Infinitely far ahead of the planing body, the flow is uniform 
and parallel to the undisturbed free surface, while infinitely far 
aft, a constant downwash, у.о is achieved. The downwash concept in 
airfoil theory is well established and is permissible in an ideal fluid 
without gravity. In the real fluid, frictional forces as well as grav- 
ity wave formation invalidate the discussion of a downwash infinitely 
far behind; however, what happens physically behind the airfoil or 
planing surface is not important as long as the mathematical formula- 
tion of the phenomena occurring at the body is mechanically compatible 
with the assumptions made about the fluid. 

With the established concept of a downwash infinitely far aft of 
the planing body, and no downwash infinitely far ahead, it can be 
seen that there will be à change of momentum in the negative Z direc- 
tion. Momentum flux resulting from the forward spray occurs at the 
control surface cutting through the spray, c d on Fig. 2. The time rate 
of totalmomentum change gives a force in the Z direction, and this 
force on the fluid is applied by the planing body. 

Let us now consider the rectangular flat surface of finite beam 
planing at velocity V. If we fix axes in the planing body, the flow 
picture is that shown in Fig. 2. The flow about the plate can be re- 


solved into two components, viz., a component along the plate and one 





transverse to the plate. It is the transverse component of the flow 
which for the moment becomes important. The transverse component of 
flow gives rise to a virtual mass described by Wagner to be the mass 
of a semi-circular cylinder of fluid. The virtual mass is not a 
physically existing mass of fluid but it does take on physical signi- 
ficance. Appendix I describes how the virtual mass concept comes in- 
to existence, but for the present, suffice it to say that Wagner has 
found the virtual mass associated with the transverse component of 
flow about the planing body to be one-half that associated with the 
lifting line. M. M, Munk, Ref. 15, shows how one finds the virtual 
mass of the lifting line by association with the potential flow solu- 
tion of transverse flow about a line immersed in the fluid. Munk's 
work is further discussed in Appendix I. Extension of the transverse 
flow about a line to the case of the planing plate is almost trivial, 
and the results lead to the mass of a semi-circular cylinder of 
diameter equal to the span. 

The force of the plate on the fluid results from the time rate 
of change of the product of the virtual mass and an appropriate velo- 
city. Similar to Weinig's low aspect ratio airfoil theory, the virtual 
mass associated with large values of wetted length-beam ratio planing 
surfaces becomes the mass of the semi-circular cylinder plus the mass 
of fluid found by Lsint where £ is wetted length and c the angle of 
Emm (see Fig. 1). 

Justification for using the virtual mass similar to that used by 


Weinig, and justification for describing the virtual mass as Weinig 





has done, may be found in Appendix I. 
From the above analysis, we find the time rate of change of 
momentum from the inlet face to outlet face of the control surface 
to be: 
Р, = mh А Б vertical component of spray momentum (2) 
where | eae | 
ћ EST +46 sinc)0V 
One should note that the Р, in equation (2) is the force of the plan- 
ing body directed vertically down (in the negative Z direction). То 
find the actual lift on the planing body, Fos it is necessary to 
eliminate the spray force. Let us define м. as the mass rate of 
flow in the spray. Then 
mm, = bS 
where § is the spray thickness (see Fig. 2). 
The lift equation by application of the momentum theorem then 
becomes: 


2 
L= (B+ gb sine) evv, - €1bg Y sine (3) 


If equation (3) is divided by 2V we have 


2 V. 2 š 
L Я те V b E 
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Transposing the spray term of equation (lh) to the left side, then 
dividing by the coefficient of v, „42+ we have: 
2 у. 
L - ev b in m 
.b$si Tw io (5) 
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Defining the length beam ratio as: 
12 
p A 


and the lift coefficient, С. аз: 


L 


equation (5) becomes: 


Cr. + 2 54 pes ie 
(A + lı sind К 
It now becomes necessary to discuss у 5.121. At a great distance ahead 
of the planing body, we found the flow parallel to the free surface or, 
in other words, the downwash velocity zero AE, 0). At the center 
of lift of the planing body, we assume the downwash to be some frac- 
tional part of the final downwash velocity infinitely far aft of the 
body. For the moment, let us call the downwash velocity at the body 
k Voz where k is a fraction. From the momentum equation (neglecting 
spray), the force on the body is: 
Е = -№ u 
and the work per unit time done by the body on the fluid 
Work = F k У; со 
The kinetic energy per unit time associated with the downwash in- 
finitely far aft is: 
2 


У. 
і оо 


КЕ = 


No] Be 


while at infinite distance ahead it is zero. The change in kinetic 


energy isequivalent to the work done so that: 


AKE = Š v,“ = ñ k v. ч 
40 +0 
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From this analysis, it is evident that k must equal 1/2 so that at 
the center of lift of the planing body the downwash is one-half its 
value at an infinite distance aft. Now it should be clear that v. /2V 
is the slope of the flow at the center of lift of the planing surface, 
and the angle whose tangent is v, 42V is called the induced angle of 
attack or induced angle of trim. Symbolically: 


tang. = m 


Equation (6) now becomes: 

Cr, 2 2 sine 

A tol (7) 

(Г + l| sin 7) б 
An inspection of Fig. 1 shows the relationships between the induced 
angle, 0. > the geometric trim angle, Y, and the ambient flow angle, Ls 

== = x + e. (8) 
Let us now return attention to equation (1) which is the two- 
dimensional solution of the planing flat plate as developed by 
H. Wagner. The angle of trim % in the two-dimensional case is 
equivalent to the ambient flow angle, e s since the induced angle 
is zero when gravity is not considered. Replacing % by &, in equa- 
tion (1) and assuming oL sufficiently small, the following approxima- 
tions can be made: 
cos ol “1 


1-cosu 


—— M is fractional and its natural logarithm is negative. 
2 созо, 


1—-cosx 
- (1-c0sog Jn aos 13 positive but very small. 
e 





12, 


The denominator of (1) then becomes: 
1 + cos Z 2 
and the lift coefficient becomes: 


E 2 sin oL (9) 
2 + sino, 


Ст, 


Now it becomes necessary to perform some algebraic manipulation 


upon equation (9). 
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o = arc sin 72-0) 


Expanding the arc sin in series: 


20, 8с, 2 
aa omm л ee 
e T2-G) B. Cu 


For values of C: « „85, the cubic term is less than 3% of-the first 
term and subsequent terms are much smaller. Neglecting the cubic and 


all subsequent terms of the series we have: 


EC 
Ric oo ENS c (10) 
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From Ref. l an expression for the spray thickness is given as 


follows: 
l r 1+ cosy _ Za l- cost, xsin Y (11) 
1 - cosg 2 cosU 1 - со 


Referring to Fig. 2, one notes that r is measured from the trailing 
edge to the stagnation point while jJ is measured from the trailing 
edge to a line tangent to the spray root. The difference between r 
and £ is considered insignificant. If in equation (7) we now let 


tana, zs w;, equation (8) becomes 








2 sin € (12) 
2 0 +i - On l- cosy , sin 
UL t x = I- cosy 2 соѕ т - COSY 
m2 - 79, (7/2X *l sin?) 
The following substitutions are made to simplify the algebra: 
Let ne =A 
1 + cos 1 = созу, w sin g Т 
T- cosv £! cost” 1 - cos? 
and GE + |1 $1197 = 
Equation (12) then takes the form: 
En 2C, | Сү + А 
Clearing fractions and collecting like powers of Ct we find: 
-= (BY+ 5B + 2 - A)cy - 2(A - BY) = (13) 
and solution by the quadratic formula yields: 
В(2+ 5) +2 - А B%+ Š B + 2 - A|° 
en d 5 + 2(4 - BY (14) 





То facilitate numerical calculations, one notes that A can be simpli- 


fied considerably. Ref. 1 contains a numerically evaluated plot of 





1h. 


$/r vs. T. From zero to about lib degrees, the curvature is slight so 
that a straight line approximation will not introduce material error 
insubsequent computations. The linearization is carried out by making 
the straight line coincide with the curve at zero and 25 degrees. The 
slope Kofthe straight line is then found per radian to be approximately 


0.2 so that: 


usos approximately 
A * (2) (0.29? = 0.1 (15) 

Values of À and B are tabulated in Table I,then the quadratic equation 
is solved. The values found from solution of equation (1h) without 
linearizing §/2 and those found when equation (15) was substituted 
for A show no appreciable difference; therefore, the values of A in 
Table I are those found using equation (15). Table II contains the 
numerical values of equation (1h), and Fig. 3 is a plot of these values, 

Modification of the two-dimensional flat-plate-airfoil solution 
for the low aspect ratio case by the method used here led to results 
which, compared with experimental data, overestimated the lift. Simi- 
larly, the planing results developed here overestimate experimental 


results. 


Application of the Momentum Theorem -~ Method II 
The fact that planing lift results are overestimated by the above 


development leads one rather logically to the conclusion that the vir- 
tual mass chosen was too large. Reflection upon the phenomena occur- 


ring at the plate then leads one to conclude that a new definition of 





15. 


virtual mass based upon the ambient flow angle rather than the geo- 
metric angle of trim should be used. 
Redefining the virtual mass, we obtain 
L = ex E + sinOl )v,__- evóp $ sin % (16) 
e' i e 


Dividing by 2V, transposing, and simplifying as before, there results: 


aa 
Vico Cr, + 2 sin % 
(sk + sine) 


Substituting for sin oc, the value found from the two-dimensional 


solution, we find: 
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Again substituting in equation (8) there is obtained: 
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Simplification of equation (18) using the linearization for §/¢ as 


previously described leads to the cubic equation: 
3, [8 ax 1,16 2 
o + | EWE nr 


[r pr Ea Jo, - Fé o (19) 


Substitution of X^ = O in equation (19) recovers the two-dimensional 


results. 
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Solutions for equation (19) are tabulated in Table III and 
plotted on Fig. |. 

The values of Cr obtained by this method underpredict the ex- 
perimental values above values of X of about 1.5 to 2. Опе may con- 
clude then that the virtual mass selected in this manner is too small 
for the low aspect ratio case and too high for the large aspect ratio 
case. The low aspect ratio airfoil treatment by this method also 


underestimates the lift. 


EXPERIMENTAL INVESTIGATION 


Purpos e 


A comparison of experimental data with the results computed by 
the two formulas for the lift coefficient developed herein was under- 
taken. It soon became evident that existing published data did not 
cover the entire aspect ratio-angle of trim range desired; therefore, 
the necessary tests were performed in the towing tank to supply the 
otherwise unavilable data. То this end, the facilities of the sea- 
plane tank at the Experimental Towing Tank, Stevens Institute of 


Technology were employed. 


Equipment and Procedure 


The model used for the tests was machined from half-hard brass 
bar stock, and provision for supporting the model in the towing car- 
riage was made by mechanical linkages shown in Fig. 5. The overall 


length of the model was 18 inches, the beam 2 inches, giving the 





174 


general appearance of a rectangular ski. The lengthwise edges of the 
upper surface of the ski were beveled at an angle of 30 degrees as 
shown in the sketch of Fig. 5. The upper surface bevel provided 
sharp edges for a clean break-away of flow. The flat bottom which 
was the wetted surface during the tests was polished and marked with 
Dykem Steel Bluing to permit ease of reading the wetted lengths from 
the underwater photographs made during each test run. The marking de- 
tails are also shown in Fig. 5. To eliminate the adverse effects of 
spray on the carriage, it was found that vertical separation between 
the model and carriage supporting linkages was necessary. To provide 
the required vertical separation, two mounting rods, approximately 
3-1/h inches long were tapped into the upper surface of the model 
near the after end. 

The seaplane tank where the tests were conducted is 313 feet 
long, 12 feet wide, and 6 feet deep when the water level is just even 
with the opening in the standpipe inside the tank. A maximum carriage 
speed of 50 feet per second can be attained with the selected speed, 
amplidyne controlled, carriage driving motor. A tachometer generator, 
the armature of which is driven by the carriage driving motor shaft, 
supplies electrical information to the field of the amplidyne system 
to provide speed control. With a selected voltage, representing a 
certain speed, the carriage speed can be maintained to within 0.06 
feet per second. A General Electric photoelectric tachometer records 


the carriage speed. 
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The camera for making underwater photographs of the wetted sur- 
face was a specially adapted Beattie Varitron, Model E, Data Record- 
ing Camera. The Varitron is a low speed, pulse type camera, the opera- 
tion of which is completely automatic. An electrical signal operates 
the shutter of the camera, and as the shutter opens, it triggers the 
photoflood lights to make the exposure. Following the exposure, the 
shutter is automatically disconnected while another circuit provides 
power for winding the film in the magazine to the next frame. When 
film winding is completed, the shutter circuit is re-engaged, and the 
apparatus is ready to make another exposure. A recording chamber on 
the side of the camera automatically records sequence numbers and in- 
formation on a platten to identify the test. The camera for the tests 
described herein was housed in a watertight box and pressurized to in- 
Sure against leakage. A control for remote focusing after the camera 
had been positioned in the tank was also provided. A light source, 
mirror, photocell, and thyratron were used to actuate a relay which 
delivered the electrical pulse to the camera shutter. 

The model was installed in the overhead towing carriage as shown 
in Fig. 5 with provisions made to set the angle of trim as desired. 
With the trim angle set at zero degrees, adjustments were made to in- 
sure that there was no angle of roll or yaw. By gently lowering the 
ski into the water and observing the bottom surface with a mirror held 
underwater below the ski, a simultaneous wetting of the entire bottom 
surface indicated proper alignment in roll. Yaw was checked by plac- 


ing a bar with reference marks inscribed thereon at the front edge of 
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the ski, perpendicular to the direction of motion and in the plane of 
the bottom surface of the ski. Then, if by moving the carriage for- 
ward, the edge of the ski remained on the reference mark, yaw angle 
was considered zero. A spirit level protractor was used to set the 
angle of trim. 

With the desired angle of trim set, a calibration in heave was 
made so that during the runs, a first order approximation on the 
wetted lengths could be ascertained. To accomplish the heave cali- 
bration, the ski was again gently lowered into the water and observa- 
tions in the underwater mirror were made similar to those for the roll 
adjustment. When the trailing edge of the ski just touched the water, 
the heave reading was taken on the large scale as shown in Fig. 5. 
This reading was listed on the data sheet (Table IV of this report) 
as the heave reference, and this reference value was subtracted from 
the heave reading made during the run. The difference (or sum of ab- 
solute values in the case of a negative heave reference) gave the 
vertical projection of the wetted length of ski excluding the wave 
rise. Then by an empirical wave rise factor, described in Ref. 18, 
the approximate total wetted length could be found. Let it be empha- 
sized here that the heave readings were used in the tests as a control 
measure only. Since there are wave disturbances in the tank after the 
first run, slight transient effects limit the accuracy of heave read- 
ings. The procedure for obtaining approximate wetted lengths from 


heave readings was as follows (see Fig. 6): 
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y 
Let: = = wetted length/beam ratio of the ski neglecting 
wave rise 
h = difference in heave (running heave minus reference 


heave, inches) 
= empirical wave rise factor = .3 


= trim angle (degrees) 


> 
H 

$ 
li 


total wetted length/beam ratio 


А л у h 
Ab? b sinT A 
E 3 h = (À - .3Xb sin T) 


Anticipating specified wetted lengths made it is possible to 
determine approximately the desired value of h so that during the 
runs, the values of h served as a check on actual wetted lengths 
being obtained. Adjustment of the load to get the desired value of 
h and, thus, the actual wetted length, was then practical. 

In order to eliminate all effects of gravity on the lift co- 
efficient obtained by the tests, all runs were made with a speed co- 
efficient, Сү» at a nominal value of twelve. Theoretical justifica- 
tion for independence of gravitational effects on the total lift coeff- 
icient at high Froude Numbers will be treated later. (Refer 
to Fig. 7). Since the beam of the model was two inches, a carriage 
speed of 27.8 feet per second was found to be appropriate. (Note that 
Су = VAB or V = Cy VED). 

For each run, a pre-determined load was applied to the model 


which was allowed to run free in heave with a fixed angle of trim and 
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at constant speed. The wetted length for each run was photographed 
so that correlation with load could be made for data reduction. 

Table IV, which is a copy of the laboratory data sheet, contains 
all of the information pertinent to each of the test runs. Table V, 
containing wave rise "correction" factors computed from the heave 
readings and the underwater photographs, was made as a check on the 
previous empirical wave rise factor and is intended to supplement 
present existing information for future tests which other investiga- 
tors may perform in the aspect-ratio-angle of trim range covered 


herein. 


Experimental Results 


Reduction of experimental data from the photographs taken and 
carriage speeds recorded are found in Table VI. 

Since most of the experimental data collected from the tests 
here were beyond the range of existing data, it was not possible to 
evaluate the reproducibility of such results, particularly in the 
limited time available. The general trend of curves faired through 
the test points obtained in this investigation is to lie between the 
curves of Fig. 3 and ) which leads one to believe that the results are 
reasonable. It is noted from the ratios of computed values based on of, 
to experimental values tabulated in Table VII and shown on Fig. 8 that 
the data collapses fairly well for trim angles of 18 degrees or less. 
Above 18° the dependence of experimental values of the lift coeffi- 


cient upon the angle of trim manifests itself, but it is not clear 
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from a physical viewpoint just what the relationship between A, TY, 


and Cr. should be. 


COMPARISON OF THEORETICAL RESULTS WIT TAL DATA 


General Discussion 

In order to determine the validity of calculations made by cor- 
recting the two-dimensional planing solution for the effect of finite 
Span, a comparison was made with experimental data. The principal 
source of data used was from Weinstein and Kapryan, Ref. 19. Since 
the existing data available did not cover the complete angle of trim, 
wetted length-beam ratio range desired, additional data were obtained 
experimentally and are presented herein. Of the several methods avail- 
able for comparing data, it was decided that ratios of computed values 
to experimental values plotted against X with vas parameter would 
give the most enlightening results. Values of the ratios using both 
method l and method 2 previously described for computing values of 
Ct. are found in Table VII and are plotted on Fig. 8 and 9. Since the 
ratios are not of the order of unity, it becomes necessary to make an 


empirical correction to the computed values in order to predict the 


lift coefficient satisfactorily. 


Empirical Correction to the Curves of Cy Computed by the Cubic Equation 


For m defined as a function of х, it is seen that as À becomes 
greater than about 2, computed values of Cr from equation (19) under- 


estimate the experimental values. The ratios of computed Or, to ex- 
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perimental C, values found in Table VIII and plotted on Fig. 8 show 
that for trim angles of 18 degrees or less, the ratio plot tends to 
collapse into a single curve. For the case of Cr, computed on the basis 
of equation (1l) the plot of ratios, Fig. 9, does not collapse as neat- 
ly as the plot of Fig. 8. For this reason, empirical correction to the 
curves of C, computed by equation (19) was elected. 

The curve shown on Fig. 8 was faired to fit most of the points 


and the empirical factor found as follows: 


E computed 


С Е 
L — - 1 (20) 


re n = empirical correction factor. 

For и. Ү = f(A) and we call this factor ny - For 
189 < # <30°, N = £(A,T) and this factor is denoted by No 

There are, of course, several suitable functions to fulfill the 


requirements on E but the ones chosen here are: 


T NN 
Nz, = 1.359 - ба (=—р—) (21) 
Mo + 125 tanh AS (22) 
From Fig. 8, the intercept value of the curve 15 recovered from YN 
when À = O, 
9 _ 180 
In equation (22), Aris defined as Y TE radians. Then when 
N = و0‎ 1, = x: 


Ihe empirical factors, Ny and T have been applied to values of 


Cr computed by equation (19) and the results tabulated in Table VIII. 


Fig. lO is a plot of the empirically corrected values ої Сү апа is 
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the one to be used for estimating the lift coefficient when gravity 
is not considered. 
Fig. ll shows the empirically corrected curves of Cy VS. À with 


experimental data points plotted thereon. 


IHE EFFECT OF GRAVITY 





One can recognize that the influence of gravity in the case 
of a planing body can manifest itself in two separate ways. One 
part of the gravity contribution is a static force which is the 
Archimedian lift or buoyancy; the other is a dynamic effect which 
attends the phenomena of the development of waves. 

Until recently, the effects of gravity on the dynamic lift have 
been practically ignored in planing lift calculations, partially be- 
cause the mechanism of the dynamic effects has not been clearly under- 
stood (Ref. 5, 17, 20), and partially because the theory, including 
the dynamic effects, is quite complex and difficult to evaluate. 

In 1937, L. I. Sedov, Ref. 2, presented a two-dimensional asymp- 
totic solution to the planing problem, including the effects of 
gravity on the dynamic lift in the region of high Froude numbers. 

For more practical applications, Hajime Maruo, Ref. 3, published in 
1951 a two-dimensional solution to the planing problem, and followed 
this in 1953, Ref. 21, with a solution for the finite span. 

it is considered appropriate to summarize here the Maruo paper, 
Ref. 21, which reveals the effects of gravity on the lift of the 


planing flat plate of finite span. 
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Maruo assumes a form for the velocity potential function Ø which 
involves an unknown — distribution over the planing surface. 
e Бе -Ю+в(% 3) tg 
p al) de dß АР. Ста ат 
0 = kK ixu * uA 
Ihis potential function is — to satisfy the free surface 
condition of constant pressure except on the planing surface. Then, 
from the potential function, a general formula for the downwash is 


found. 


karik ofa ihe Dee e Denon Done} 


al [| ра е En torpe El 


Maruo points out i the downwash can be separated into two com- 
ponents; one corresponds to a profile characteristic and is present 
even in the two-dimensional flow; the other is a downwash arising 
from the effect of finite span. 

The downwash condition applied at the planing surface to satisfy 
the condition of no flow through the plate, gives rise to an integral 
equation for the unknown pressure distribution over the plate. 

Maruo, by a procedure analogous to Glauert's solution of the 
finite span airfoil problem, simplifies and solves the integral equa- 
tion for values of \ <1/2, or for aspect ratios greater than about 
two. Fig. 7 is a graphical representation of the solution when A-—>0, 

One should note on Fig. 7 that there are two curves, I and II. 
Curve I is the plot of a numerical solution to the exact equation 
while II is a plot of the simplified equaticn evaluated analytically. 


Attention is directed to the shape of the curve of Fig. 7 which shows 





that the theoretical value of Cy not considering buoyancy effects, 
rises from a Cy of 0.75 and approaches asymptotically a constant value 


at higher C Most plots of experimental data show the total lift co- 


у: 
efficient without subtracting the coefficient of buoyancy 8/2 €V^b | 
from the data; therefore, the effect of gravity on the dynamic lift 
appears to be a monotonically decreasing function of Cy rather than 
an increasing function as shown by the theory. 

The practical results of Maruo's work are given by Equation (h2) 


of Ref. 21 which, when transcribed to the notation used herein, be- 


comes: 
wee 

1 

M+ = 

y 3 





ot, 


where y is a function of the speed coefficient, Gr. 
A graph of the y/function, given in Ref. 21, can be linearized 


to the following approximation: 


em 
y - 0.667 + ur 


Then Cr becomes : 


2 
ООШОТ 07333 + 0.233 (23) 
C 2 
+ 


From equation (23), it is evident that the lift coefficient becomes 


independent of the dynamic effects of gravity for 


0. TX 
т. 


One should be aware of the fact that the mathematical validity of 


equation (23) was based on \ <1/2. 
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For small X (such that A<< 1), it is seen that the first term 
of the binomial expansion of the inequality yields C, O. 6VÀ which 
agrees quite well with the empirical value of Ci /VA > 1 established 
by Paniel Savitsky, Ref. 18, as a criteria for independence of gravity 
in the flow. 

The effect of gravity on the dynamic lift from this analysis is 
seen to become negligible at low values of Cy. Those investigators 
who insist that the value of Су above perhaps lO is necessary to elim- 
inate the gravity effect refer to the combined buoyant and dynamic 
effects. The dynamic effects become small at low values of Сү» but 
until the dynamic pressure, 1/2 ет", becomes large enough to make the 
ratio of static buoyant force to dynamic pressure a small percentage 
of the total lift, buoyant effects will persist. 

For most practical cases, the lift coefficient, as Су approaches 
the order of 10, is clearly independent of all gravity effects. If 
one subtracts from the total lift coefficient the buoyant 1ift con- 
tribution defined by PR it is seen that the remainder is in- 


dependent of all effects of gravity only for Ci > 0.6VA when А < 1/2. 


SUMMARY AND CONCLUSIONS 
The results of this investigation show that theoretical evalua- 
tion of the effect of a finite span on the planing surface problem is 
quite complex. To remove some of the complexity, the problem was con- 
sidered in two phases. First, the effects of gravity were neglected, 
and an ideal fluid assumed; then second, the effects of gravity were 


discussed. 
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Neglecting gravity, equation (1h) was obtained from a combination 
of the two-dimensional solution and application of the momentum 
theorem to a selected virtual mass. Equation (19) was found in a 
similar manner but with a different virtual mass. Both equations (1) 
and (19) overestimate the lift coefficient for aspect ratios greater 
than about one-half which indicates that the choice of virtual mass 
for that range of aspect ratios is too high. Equation (19) underesti- 
mates the lift when wetted length-beam ratios are greater than about 
two (aspect ratio less than one-half), while equation (1h) overesti- 
mates the lift for all length-beam ratios. 

Fig. l shows cross-sectional areas of the "volumes" containing 
the virtual mass for each case. Tables II and III give tabulated 
values of the calculations made by each of the above methods, and 
Table VIII gives the empirical corrections to make values computed 
by equation (19) agree with experimental data.  Empirically corrected 
values of the lift coefficient are plotted on Fig. 10, and comparison 
of these curves with experimental data are shown on Fig. 11. 

While the analytical results obtained are considered reasonable, 
empirical correction to the theoretical analysis was required to get 
better agreement with experimental data. Thus, the development can be 
considered as a semi-empirical method for evaluating tne effect of a 
finite span on the solution to the planing problem in a weightless 
O. 

For the analysis made independent of gravity at length-beam- 


ratios greater than about 2, the method leading to equation (19) 
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indicates that the virtual mass may be too small. The virtual mass 
associated with equation (ll) may be too large. The most prominent 
physical explanation for the inconsistency between theoretical and 
experimental results is that the downwash along the wetted length of 
the body is not of constant magnitude. Bollay, Ref. 16, indicates 
that free vortex lines are shed along the chord of the low aspect 
ratio airfoil so that at the trailing edge,the contribution to down- 
wash there is not only from the vortex wake, but also from the shed 
vorticity ahead. This, in effect, produces a curvature of the flow 
from leading to trailing edge which has the effect of introducing an 
apparent camber to the foil. 

For the low aspect ratio (large A) planing surface, the curvature 
effect may have a significant influence on the lift, and is offered as 
a possible explanation for the difference between predicted and experi- 
mental values. 

lt is worthy of notice here that the technique similar to that 
used for modifying the infinite span airfoil for the effect of finite 
span is a logical procedure to use for the planing surface. As 
pointed out by many other investigators, the flow for most practical 
purposes can be considered ideal and weighless in the planing regime, 
and from the stagnation point to trailing edge on the lower surface, 
both airfoil and planing body have almost identical pressure distrib- 
utions. 

The consideration of gravity led to results which show that the 


effect of gravity on dynamic lift is a function of the speed coefficient 
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and the wetted lengtn-beam ratio. It was found that for all practical 
purposes, the effect of gravity on dynamic lift could be disregarded 
ШКО, > 0.6VÀ when À—21/2. Buoyant effects remain significant 
until the dynamic pressure reaches such a value as to make the co- 
efficient of buoyancy very small compared with the dynamic lift 


coefficient. 
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APPENDIX I 


General 
Ihe use of incompressible potential flow theory to treat the 

effects of low aspect ratio on wing lift has been carried out for 

the case of a flat plate airfoil in at least two distinctly different 

ways, but with essentially the same result. The two different methods 

cited are: 

(1) The solution by F. Weinig, Ref. 1h, whereby the momentum theorem 
is applied with an estimate of virtual mass, and an estimate for 
the induced angle of абфаск, о ., is made based on cascade theory. 

(2) The integral equation solution of W. Bollay, Ref. 16, in which 
the boundary condition—that there be no flow through the 
plate--is satisfied by equating the transverse velocity com- 
ponent of the stream, V sing, to an integral expression estab- 
lished for the downwash at the wing. 

It should be noted that both of these methods are developed for 
an airfoil in an incompressible fluid neglecting gravity; therefore, 
it is reasonable to expect the results to be analogous to that рог- 
tion of the planing problem which is treated independent of gravity. 
Wagner, Perlemuter, Perry, and Shuford, Ref. 1, 5, 6, 12, 20, and 23, 
have indicated that the integral of pressures over the wetted surface 
of the planing body, excluding the region about the leading edge 
forward of the stagnation point, is almost identical with the pres- 


sure distribution over the same portion of a flat plate airfoil; 
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thus, the planing surface-airfoil analogy seems to be reasonably 
founded. 

It is not the purpose of this discussion to present the details 
of the airfoil solutions of Weinig and Bollay, but rather to show 
how they lead to a better understanding of the virtual mass chosen 
to modify the two-dimensional planing surface solution for the effect 
of a finite span. The fact that the solutions of both Weinig and 
Bollay lead to essentially the same result is graphically revealed 
in Fig. 12. The experimental points shown in Fig. 12 were obtained 
by H. Winter, Ref. 2l, and are included to show the relative merit 
of each of the theoretical results displayed. One obvious reason for 
the difference between the theoretical solutions of Weinig and Bollay 
is that mathematical complexity forced Bollay to make several simpli- 
fying assumptions in order to solve the integral equation established 
in his formulation of the problem. Although Bollay's treatment of 
this problem seems to be established on a more rational foundation, 
Weinig's results correlate better with experimental data as pointed 


out by H.R. Lawrence, Ref. 25. 


Virtual Mass 

Modification of the two-dimensional planing surface problem for 
the effect of a finite span by application of the momentum theorem 
makes necessary the choice of an appropriate virtual mass. Attention 
is directed to the low aspect ratio airfoil solution by Weinig wherein 


the virtual mass selected is composed of the virtual mass associated 
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with a lifting line plus the streamwise projection of the chord 

bl sing. In order to better understand Weinig's choice of virtual 
mass, a review of the salient features of the virtual mass concept 
as given by M.M, Munk, Ref. 15, are considered. A brief review of 
the most important steps leading to the definition of virtual mass 
follows. 

If an incompressible fluid at rest is given an impulse which re- 
sults from the action of normal pressures alone (frictional forces 
assumed zero) the resulting flow is potential flow provided no dis- 
continuities exist and all spacewise partial derivatives of the poten- 
tial function are continuous. If it is assumed that the impulsive 
force creating the flow acted for some finite distance, work would be 
done on the fluid, and since the fluid was originally at rest with 
zero kinetic energy, the work done represents the kinetic energy 
stored in the flow. It is convenient to replace the total energy 
of the fluid particles moving in many different directions with the 
valent energy of a certain geometrically defined mass of the fluid 
moving with the velocity of the body. This assumed mass, introduced 
for mathematical convenience, is labeled virtual mass and may be 
different in magnitude for different directions, depending on the 
velocity chosen in its definition. 

If avery thin airfoil with axes fixed in the body moves with 
a steady velocity in air which was originally at rest, velocity about 
the airfoil at angle of attack,& , can be resolved into two components, 


one normal and one parallel to the foil as shown in the sketch. 
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In an ideal fluid (no friction), 







the component parallel to the air- 
EE foil does not contribute to the 
= forces. Attention is then concen- 
q trated on the normal component which 
is referred to as the transverse 
flow component, The line pq on the sketch represents an element of 
the airfoil of span b. (It has been shown that when the span to chord 
ratio is large, the entire foil may be represented by such a line.) 
The wake at some distance behind this airfoil behaves like a flat plate 
of span b moving with velocity v perpendicular to the plane described 
by the plate. The complex potential for the flow net about this hypo- 
thetical flat plate is given by: 
w= + iy = iv (z - №“ = (50°) (А1-1) 
where 2 = х + iy, у = transverse velocity component, and b/2 = semi- 
span. The product, P 6, can be thought of as the impulsive pressure 
necessary to create the flow field in which kinetic energy is stored. 


Then, the expression for kinetic energy in terms of the potential func- 


tion is: 
d 
T = ме] в 0 ag ds (А1-2) 
where © = fluid density 
dg/dn = component of transverse flow velocity 
ds = surface element over which the integral is performed. 


This is the total kinetic energy of the flow, and as previously 
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stated, it can be expressed in terms of the virtual mass, m, as 
Т = 1/2 mv (A1-3) 
where v = d@/dn = constant on the surface of the plate. The total 
kinetic energy of the flow as found by the potential function is 
equated to the energy expressed by the virtual mass so that 


t ея - 
Š 


mv? (Al-Iı) 


NIH 


Since pus constant and - v = constant, there is obtained: 
late 
vol = 4, | g Was (A1-5) 
— 


This defines the volume of fluid of density p associated with the 
transverse velocity. It should be noted that the virtual mass is 
not à physical body of particles; its properties, as stated previously, 
are directional such that for a direction other than that of the trans- 
verse flow, the virtual mass may be entirely different. 
On the surface of the hypothetical plate (i.e. y = O in the com- 
plex variable notation): 
2 
ERR х2 (A1-6) 
Substituting the value of Ø from equation (A1-6) into equation (A1-5), 
and noting that або is constant over the plate equal to v, the 
volume of fluid defining the virtual mass is obtained: 
Vol 3" > -xf dx (A1-7) 
Since the surface over EON S integral is performed extends 


from -b/2 to +b/2 on both upper and lower surfaces of the plate, the 


final expression for the volume becomes 
b/2 s p? 
Vol = 2 2 -x ах = (A1-8) 


-b/2 
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This is the volume of fluid containing the virtual mass and is shown 
here to correspond to the mass of a circular cylinder of diameter 


equal to the span of the plate. 


The Lifting Line 

In the previous derivation of lift over the span b and an element 
of length dl, the downwash velocity, v = d@/dn was necessarily con- 
stant along the span and was assumed to have no chordwise variation. 
The same constant downwash velocity is known to exist in the case of 
the vortex system consisting of bound and trailing vortices such that 
the circulation of the bound vortex varies elliptically along the 
span. The bound vortex is known as the Lifting Line, and it can be 
taken to represent the entire wing when the ratio of span to chord 
is large. 

For multiple lifting lines disposed vertically in close proximity 
to each other, Munk has indicated that the interference existing be- 
tween the flow patterns associated with each line reduces the virtual 
mass of the single line. The mathematical complexity involved for 
the case of three or more lines limits the discussion here to only 
qualitative, order of magnitude results. 

In the case of a biplane, two lifting lines separated by a small 
vertical distance, h, can be assumed to replace the wings. The re- 
sulting volume containing the virtual mass associated with the trans- 
verse flow is given by Munk as approximately the cylinder of diameter 


equal to the span plus the rectangular mass enclosed between the upper 
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and lower wings: 2 
Vol d + bh 
The Effect of a Long Chord 

The virtual mass associated with the lifting line, used in 
application of the momentum theorem to find the lift for a finite 
span airfoil leads to the prediction of reasonable values for span 
to chord ratios which are large. As the chord lengthens with respect 
to the span, the single lifting line representation becomes invalid, 
principally due to Be that the free vortex lines no longer 
leave the trailing edge in a flat plane of the main stream, but in- 
stead are shed along the chord forward of the trailing edge to form 
a three-dimensional vortex array in the wake. 

In order to satisfy the boundary conditions d@/dn = V since 
on the surface of the airfoil, and of the streamlines flowing smooth- 
ly off the trailing edcze, Bollay, Ref. 16,22, has assumed a distríbu- 
tion of bound vortex lines along the chord. After certain simplifi- 
cations, Bollay has solved the integral equation for the distribution 
of circulation among these vortex lines (or lifting lines). 

Interpreting Bollay's results in terms of virtual mass, the 
chordwise spacing of any pair of bound vortex lines at a distance al 
apart leads to an equivalent vertical spacing of 9@ зто. Consider- 
ing this in the light of Munk's theory of the biplane and multiple 
lifting lines, the integral along the chord of the contributions of 
all vortex pairs gives the contribution bf sina to the volume of the 


virtual mass as an approximation. The total volume associated with 
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an airfoil of small span and long chord then becomes m + bl sind. 
It is observed that this is the exact expression chosen by Weinig. 
It should be noted, however, that & in this case needs an inter- 
pretation. Weinig has assumed it to be equal to the geometric angle 
of attack of the airfoil, but in the spirit of the above derivation, 
it would appear more logical to use the ambient flow angle in the 
immediate vicinity of the airfoil rather than the geometric angle of 
flow. 

Modification of the two-dimensional airfoil solution, using the 
momentum theorem, with the virtual mass based on geometric angle of 


attack, then on the ambient flow angle, follows. 


Modification of the Two-Dimensional Airfoil Lift Coefficient for the 


Effect of Finite Span 
Method 1 - virtual mass based on geometric angle of attack,o&. 


From the momentum theorem: 
L=ñv, = net (fl + A sinoQv. (A1-9) 
ley +00 
wo? 
where m = G + b зіп) 


X d 
b 


From energy relationships: 


V. 
= = tank zw. (A1-10) 


Dividing (Al-9) by 2V, substituting the result of (A1-10) and solving 


for@., there is obtained: 


a AG, 


Б. —-—-A——— e Al-11 
= (2у) (61) = A sino) 1* UA sinog | 
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L 
щеке С. =. 
L 1 We 
s ev 
Now: 
OL =a, +O, (41-12) 
and from two-dimensional theory: 
= 2 sing, cog, (A1-13) 
Dividing (Al-13) by: 
i 
> = 2 sind, cose, = sin20& (Al-LIı) 
E 
2% ce e 7 (A1-15) 
C 
Writing the series expansion for arc sin Ma 
t 2 Eus. 


Neglecting all terms in the series higher than the cubic, and sub- 


NeuLing fran (A1-16) and (41-11) into (A1-12): 


С Cae AG 
| оао E 
= 3 me Ë 6(" 77+ I sins el 
Solving: 
3 dm s 


Values of substituted into (Al-18) show no appreciable dif- 
ference from the results obtained as follows. 
iifequation (41-13) let: 
cos x, x.l 
snol, x% 
Then (A1-13) becomes: 
= TI, (A1-19) 
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Multiplying (Al-12) by 27, then substituting from (Al-11) and (A1-19), 


there is obtained: 


eT X0, 
2 T = Cr. * T+ IX sino (41-20) 
Equation (A1-20) simplifies to: 
e MS eee (A1-21) 


L 1 + 2 IT 
п + WX sine 
Values of C, computed by (A1-21) are tabulated in Table IX and plotted 
on Fig. 13. It is observed that this method overestimates the ex- 


perimental values of Winter. 


Method 2 - virtual mass based on ambient flow angle, Of. 


Similar to method 1, there is obtained ford. : 
AG, E 
оё. = a ЕЎ (41-22) 
d EN sine. 


Approximating again: 
cos % 1 
e 
sinot «o 
e e 


and substituting o = dm the equation: 





¿TT 
OC =O, + Ol 
there is obtained: 
C AC C AC 
L L L L. 
= — + = = — + — Al- 
= 2m 7r+ loc 2m Cr, И 
W + Un 5 


Simplifying and collecting coefficients of like powers of Ст: 


s fr cube A. m— 





Equation (Al-2)) is solved by the quadratic formula, values tabulated 
in Table IX and plotted on Fig. 13 for comparison with method 1 and 
the experimental data of Winter. It is observed that Cr is underesti- 


mated by equation (A1-2)). 
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APPENDIX II 


Foundations of the Theory of Planing 


Translation of Paragraph 75 of the Russian Textbook, Water Resistance 
to Movement of Ships, by P. A. Apuchtin and J. I. Voitkunsky, 


Moscow, 1953. 


The viscosity of a fluid affects the pressure distribution along 
the planing surface very little and separation of the boundary layer 
is not observed on planing surfaces. This provides the basis for in- 
vestigation of the phenomena of planing by methods of hydrodynamics of 
a perfect fluid, in particular the potential flows. 

In the domain of the theoretical investigation of planing, the 
Soviet scientists played a leading part. In 1929, G. E. Pavlenko 
first obtained the solution for the lifting force and resistance of 
a plate of infinite span, planing at a small angle of trim on the sur- 
face of a perfect fluid. 

Later, G. E. Pavlenko developed the foundations of the theory of 
planing, taking into account the finite span, studied the physical 
phenomena connected with planing, and also investigated the question 
of recalculation of model test results of planing surfaces to full 
Size. 

In subsequent works, the two-dimensional problem of planing of 
plates on the surface of a perfect fluid was attacked along two paths: 

1) by the methods of streamline flows, not taking into account 


the weight of the fluid, 
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2) by methods of the theory of waves of small height taking into 

account the weight of the fluid. 

The investigations of planing by means of streamline flows for 
infinite depth and for finite depth of a fluid were conducted by the 
academicians S. A. Chaplygin, M. I. Gurievich, A. P. Yampolsky and 
Wagner. 

The solutions of the planing problem of a curved plate, taking 
into account the weight of the fluid were given by the academicians 
N. E, Kotchin and L. I. Sedov. 

The planing of the surfaces moving one after another (in tandem), 
interesting in connection with the action of steps, for the case of 
gravity-less fluid was investigated by L. I. Sedov. 

Let us consider the origin and the calculation of water resist- 
ance for the planing flat plate on the surface of a perfect fluid. 

Large pressure gradients are observed along the lower surface in 
the movement of the plate. Plots of the pressure distribution, 
measured experimentally along the centerline of the plate and along 


1) Along 


its edges at various angles of trim, are shown on Fig. 13 
the span of the plate pressure changes are not significant, while 
normal to the span in the region of the leading edge the pressure 
gradient is very large. A large drop of pressure ís observed at the 
lateral edges, where the pressure abruptly drops to a value corres- 
ponding to the pressure in the undisturbed fluid. 


In regions of large pressure gradient (in the regions of inter- 


section with the free water surface), spray jets are formed. The pro- 


(1) See Fig. 1) 
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jection of the reaction of the jet, in the direction opposite to the 
direction of the plate movement, represents spray resistance. 

Let us consider the flow of a parallel stream of velocity V onto 
a stationary plate. Consider the depth of the fluid limited, equal to 
h. The influence of gravity on the flow of the fluid will be neglected. 

Flow of the fluid in such a motion is steady. Applying Bernoulli's 
equation to the particles on the free surface (Fig. 135) (1) and taking 
into account the fact that pressure on the free surface is constant, 
it can be shown that the velocity along the surface of the jet is con- 
stant also, and is equal to V. 

In the actual motion of the plate, the velocity in the jet is 
double that of the plate motion. Isolate by means of a control surface 
ВАА лозеа volume of the rivia ABODEP!) Aecording to the laW of 
momenta, the flow of momenta through the control surface is equal to 
the vector of external forces applied to the volume, ABCDEF, 

The difference of momenta flow through the sections of the con- 
trol surface, AB and CD, is equal to VQ, where Q is the mass flow in 
the spray jet. 

Projecting the flow of momenta through the control surface ABCDEF 
and the vector of forces applied to the plate on the axis OX, we get 

R. w QvQ (1+соѕо) 


and the total reaction is directed normal to the plate and is equal to 





R 
уе м lecosx _ о 
S E (Pa sino ` er > 


(1) See Fig. 1l 
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These formulae are valid for any value of the angled; the theory 
discussed can be considered therefore as the non-linear theory of 
planing. 

As the fluid velocity in the spray jet in the case of a stationary 
plate is equal to V, the volumetric flow is 

0 = 15 
where Á is the thickness of the spray. 

The force, К, in the case considered is equal to the spray re- 
sistance, R6 of the plate. Finally, the formulae for computation of 
the spray resistance and total hydrodynamic reaction in planing on the 
surface of a gravity-less fluid of arbitrary length at any angle of 


trim will take the following form: 


Rç = evs (1+с0506) (136) 
В = evs ctg 5 (137) 


If the plate of infinite span planes on the surface of a fluid, 
located in the field of gravity, waves are formed in its wake. In 
this case, the additional force of wave resistance will act on the 
plate. 

In the case of the motion of a planing plate on the surface of 
a perfect fluid, therefore, the total resistance is the sum of spray- 
making resistance, Re, and of the wave resistance, Rn. Considering 
formula (135) (1) in the light of the above derivation, we find that 


Ry * Rc - D tan«, 


(1) R, = D tanol+ Rip (135) 
D = total lift, i.e. dynamic + Archimedian 
R = the force of friction 


tp 
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The sum of forces Rp + Re represents the projection of the re- 
sultant force of normal pressures, distributed along the plate, onto 
the direction of motion, i.e. represents resistance due to pressures. 

With the growth of relative velocity, the part played by the 
gravity forces diminishes in comparison with inertial forces; there- 
fore the part played by wavemaking resistance in the total resistance 
of the plate diminishes. 

The comparison of the relative importance of the wavemaking and 
spray resistance in the total resistance of the planing of the flat 
plate is given on Figure 136, (2) 

The complicated calculations permitted L. I. Sedov to obtain the 
formula for computation of the ordinates of the free fluid surface, 
hydrodynamic pressure, lifting force, resistance and the hydrodynamic 
moment acting on the planing surface. 

In cases of large relative velocities (small values of the para- 
meter Y = 5-1-2), the formula for calculation of the lifting force, re- 


sistance, and hydrodynamic moment of the flat plate were obtained by 


means of Sedov's theory: 


R =R =PT(1 -TY - Ва У (138) 
R. = Rx (139) 
2 
eir 8437, 2 
М = E 2 en) a. yal (1,0) 


where a = - and F = V vet. 
These formulae show that at the large relative velocities and 


small angles,& , the lifting force and hydrodynamic moment in planing 


(1) See Fig. 1l 
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are linear functions of the angle of trim. 

More accurate numerical calculations of the lifting force and 
moment by Sedov's theory were made by U. S. Chaplygin for a large 
variation іп the parameter, ə Results of these calculations are 
compared on Fig. 1371) with results computed by formulae (138) and 
(110). On Fig. 137(a) and (b) are plotted the dependence of С, fet 
and the relative distance of the center of pressure from the rear 
edge of the plate 2. / f| on the number Fr. = v/vg2. Comparison of the 
curves shows that beginning with Fr. = 2.8, i.e. withy- 0.061, the 
results of the refined numerical computations and of calculations by 
approximate formula coincide. In the case of very large Froude 
numbers, Fr., it can be assumed that ye O; this assumption means 
physically that the gravity of the fluid is not considered. With 
the introduction of such an assumption, the formulae for the flat plate 


are further simplified to: 


Е = Ё, = Ma V oc (11) 
К = опа vere (142) 
м = pra voz (143) 


Were the distance of the center of pressure from the trailing 
edge of the plate computed by means of the above formula for the 
moment, it would be found that p = (3/11)2a (1.е. Mp = 75A); this 
agrees well with test results. 

Let us compare formula (141) and (112) with the formulae for the 
computation of corresponding quantities obtained in the wing theory 


for a flat plate of infinite span in a streamline flow at small angles 


(1) See Fig. 1) 
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of attack: 
R_ = 260ITa Va (144) 
M = ¿Ia V“ (115) 
The flat plate represents the simplest wing. Тһе comparison of these 
formulae shows that in planing at large Froude numbers the lifting 
force and moment of the plate is half that obtained in the flow 
around the plate deeply submerged in a fluid. This relationship 
points to the analogy of the flow on the lower part of the plate 
in both cases. 

The development of the analogy between the planing surface and 
the wing has permitted application of a number of relationships ob- 
tained in wing theory to study the phenomena of planing. 

From the theory of the wing it is known that in a flow about a 
plate at an angle of attack0f, in accordance with the eona oi 
Joukowsky, a lifting force occurs, normal to the direction of the 
flow at infinity. The velocity at the trailing edge of the plate 
has a finite value, and the streamlines flow off smoothly. At the 
sharp leading edge the velocity tends to become infinite, and in 
connection with this, the pressure in the vicinity of the leading 
edge drops sharply.  Hydrodynamic pressures are normal to the surface 
of the plate, and their resultant, P, is also directed along the 
normal to the plate (Fig. 138,2). (1) In addition to the force Р, а 
suction force, S, acts in a direction along the plate. The total 
hydrodynamic reaction, the lifting force of Joukowsky, R, is the 


vector sum of forces P and S: 


(1) See Fig. 1h 
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E $P E 

The analogy of the flow in the space under the place in the 
case of planing at small angles of attack, and in the case of the 
plate fully immersed in a fluid, is valid from the trailing edge to 
close proximity of the leading edge. Analogous also are pressure 
distributions in this part of both flows. 

The flow is different only in the region of close proximity to 
the leading edge where the spray jet directed forward occurs along the 
plate (Fig. 138,0). 4) 

From the above reasoning, the results of the study of flow about 
a plate fully immersed in a fluid can be used for the study of flow 
under the planing surface if the flow about the leading edge (where 
velocities grow to infinity) is replaced by streamline flow reproduc- 
ing the spray jet. Such a replacement at small angles of trim will 
cause but an immaterial change in the area of the noraml pressure plot 
diagram at the lower edge of the plate. Taking into account the fact 
that the flow in planing occupies only the lower part of the space, it 
can be expected that the resultant of the normal force R at the planing 
surface is equal to: 

n sO. 
The validity of this relationship was already noted earlier. 

Developing the analogy, we can suppose that the reaction of the 
spray jet is directed in the opposite sense to the force of suction S, 
and in magnitude is equal to R, = S/2. Making use of formula (136) 


for the spray resistance in the case of planing on the surface of a 


(1) See Fig. 1) 
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gravity-less fluid 
R. = Spr (1*cosX), 
and letting coso&= 1 for small angles, and R. = S/2, we can compute 
the thickness of the spray jet: 
л 
тҮ 
where S = suction force of the wing and is given as: 
5 = 2Пра Ve 

The validity of the above analogy is confirmed experimentally at 
large Froude numbers up to values of the angle,&, from 7 to 10 de- 
grees. 

In the case of flat plate planing at an arbitrary angle of trim, 
the resultant of hydrodynamic forces, R, can be calculated by 
formula (137) 

К = evs ctg X 

The curve K = f (et), showing the ratio of the magnitude of force R, 
calculated by formula (137) to the magnitudes calculated by formula (141), 
(ive. on assumption of small angles, 0%), is plotted on Figo. 139. (1) 
Fig. 139 shows that with increase of the angle,CC, the magnitude of 
the resultant hydrodynamic pressure, computed by non-linear theory, 
becomes smaller than that computed by linear theory, i.e. the analogy 
between the planing surface and a wing is disturbed. 

The theory of planing permits one to appraise the effect of 
shallow water and of limited channel width on the hydrodynamic charac- 


teristics of planing surfaces. 


(1) See Fig. 1) 
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The detailed investigation of the effects of shallow water on 
planing upon a gravity-less fluid at various angles of trim was made 
by U. S. Chaplygin. 

As a result of these investigations it is established that, in 
planing on shallow water at small angles of attack, the analogy with 
the wing moving between parallel walls is observed. Theoretical cal- 
culations show that hydrodynamic characteristics in motion on shallow 
water depend on the ratio :h, and the angleo¢, where h is the water 
depth. 

With increase of the ratio Ü:h, i.e. with a decrease of water 
depth, the lifting force at constant velocity increases. 

On Fig. 140, 1) curves are given permitting one to estimate the 
increase of the lifting forces caused by shallowness of water at various 
values of the parameter L:h and angles of trim, oC. 

Taking into consideration the effect of channel walls on planing 
without consideration of gravity indicates that, at the ratio of 
channel width to the width of the plate b:B «7, hydrodynamic forces 
acting on the plate increase; at b:B> 7 the effect of walls leads to 
a small (less than 2%) reduction of forces as compared to motion on 


the surface of unbounded fluid. 


(1) See Fig. 1) 
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TABLE II 


Computation of ©, Using m Based on Geometric Angle of Trim, T 


— 


Ali 2 
A A Le 2-4], 2a - m 


1 2 3 7 5 6 7 

m T+ AT. +2- = 

À Be E ae 2-8 — 1/0: c 

0.2 2 3.683 13.56lı -.55T 3.606 .077 

In 3.871 1. 985 UNIT 3.722 .119 

6 li, 06lj 16.516 =e S 3.816 .218 

9 lj, 360 19,010 -2. 6l hols . 315 

12 lı. 665 A02 -3.602 li, 261 Mer 

18 5.301 285101 55633 1.713 .558 

21 5.968 35.617 -74803 5.271 69 

30 6.661 bh. 369 -10.100 5.853 ,008 

0.) 2 2.365 5.593 -3283 2.30) 061 

l 2.18), 6,170 —„583 2„363 Лоу 

6 2.608 6.802 2590) 2.129 ‚179 

9 2.801 7.816 —1 3830 2.537 . 264 

12 3.00) 9.02) -1,958 2.658 . 316 

18 3.1134, 11.792 -3.165 2.937 ‚197 

2) 3.887 15.109 -1. 500 3.256 .631 

30 li, 383 19.210 -5, 988 3.635 „718 

0.6 2 1.925 3.706 = 1.875 „050 

li 2.022 la 088 =. 100 1.920 .102 

6 2.123 1.507 — i627 1.970 .153 

9 24202 5.208 =.999 2.052 „230 

12 2.150 6.003 VAMO 3.12 . 308 

18 2,812 7.907 -20 3113 2.359 ‚153 

2), 3.210 10. 30); 3.116 2.624 „586 

30 3.628 13.133 ION? 2.918 o 706 

0.8 2 1.705 2. 907 ER 1.662 . 013 

ly 1579). 3,208 =. 309 1.703 .088 

6 1.880 3.534 0190 1.715 155 

9 2,022 11.088 ==„ 793 1„815 „207 

12 2.173 15222 -1:135 1.893 „280 

18 SEO 6.255 213052 2,079 o 422 

2l, 2.860 8.180 -2.868 2.305 „555 

30 3.21), 10.52) -3.931 2.568 .676 





TABLE II (Cont'd. ) 


1.0 


2.0 


6.0 


8.0 


2 3 L 5 6 7 
о + + 2. 
a EE e Gf 2 (A-BT) ۸ 6 +(5) En 
2 1.574 2.1177 =. 186 1.536 . 038 
h 1.652 2. 729 -25h 1.573 .079 
6 1.735 3.010 —MO7 1.612 „123 
9 1„866 3.182 ~.670 1.677 2189 
12 2,007 lı. 028 sor) 1.718 .259 
18 2.314 5.355 -1.685 1.916 . 398 
2h 2.652 7.033 -2.539 2,120 .532 
30 3.016 9.096 -3.520 2.361 .655 
2 1.310 1.716 ~ 06 1.285 . 025 
|, 1.375 1.891 — 21s 1.321 205); 
6 nnn 2.085 = „оне 1.357 .087 
9 1.555 2.118 —. 123 1.12 13 
12 1.675 2.806 ¿612 1.471 . 20h 
18 1.911 34907 1.2091 1.605 .336 
2h 2.238 5.009 — 0881 1.768 470 
30 2.561 6.559 -2.698 1.965 .596 
2 1:2178 1.388 -.036 1.163 .015 
li 1.236 1.528 -.090 1.199 . 037 
6 1,298 1.685 161 1.23) 06h 
9 1.399 1.957 -. 300 1.287 „112? 
12 1.509 2277 -.1178 1.341 .168 
18 1. 75h 3,077 — 9 1.460 .29}, 
2}, 2.031 1.125 -1.552 1.60); B n 
30 2.333 5.143 ROT 1.775 .557 
2 1.131 1.286 22027 1.122 02 
l 1.190 1.116 -,071 1.160 .030 
6 1.250 1.563 =. 133 1.196 .05) 
9 1.347 1.81) -o 259 1% 217 . 100 
12 1.453 25111 4023 1.299 „15 
18 1.692 2.863 2.562 1.11); 278 
2h 1.961 3.816 —1 42 1.550 nal 
30 2.257 5.09) -2.150 1.71.5 „512 
2 1.112 1.236 ¿027 1.102 ‚ 010 
lh 1,4107 1.362 -.062 1.110 .027 
6 1.225 1.501 =. 120 1.175 .050 
9 1.321 1.715 221230 16227, 09) 
12 1.126 2.033 _. 395 1.280 116 
18 1.661 2.759 07 1.392 .269 
2h 1,927 3.713 -1:387 1.525 e 102 
30 28219 1. 92 =2 7001 1.686 „533 
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2 


m B +2+24 


N DEM 


12 


2li 
30 


3 


1.099 
1.153 
1.211 
1.305 
1.109 
1.612 
1.906 
2.196 
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li 
1.208 
1.329 
1.166 
1.703 
1.985 
2.696 


3.633 
lı. 822 


2 (A-BT) 


=„@20 
-.057 
MEL 
22220 
-. 379 
-.796 
-1.355 
~2.0h0 


. 010 
. 025 
. 017 
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„26! 
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TABLE III 


Computation of Lift Coefficient Using a Virtual Mass, m 


+ [ - 
2 6.283 
E таз ( ` 
2 3 
a. СовГе оѓ 
L 
р 
2 -3 „354 
4 T1 
6 -% 204 
9 -3,266 
12 -5 228 
18 -3.1 
24 -3 077 
30 -3 002 
2 -2 ,826 
4 22 774 
6 ara 
9 22612 
12 -2.563 
18 -2 „406 
24 л 
30 -2 .090 
2 -2.563 
4 -2 „496 
6 -2 „430 
9 Ec 
12 EDO 
18 2052 
24 -1.833 
30 -1.635 
2 -2 „475 
4 -2 „404 
6 -2,333 
9 2209277 
12 E 1 
18 -1.908 
24 -1.696 
30 -1,484 
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Coeff of 


4 


Cy 


q 


6.059 
Ое 
62178 
6967 
6.556 
6,554 
6.712 
62589 


4.950 
4.900 
4.850 
4,774 
4 „698 
4 ,547 
34 597 
4,245 


4,395 
4,290 
4 „185 
4 ¿027 
3.870 
5,099 
5.258 
22928 


4,210 
4.086 
5.969 
d 
3.092 
5 220 
2 „852 
2,481 


Based on Effective Angle of Attack, 0. 


1 2 
) L = 2 679 x C, 
- же с. - 


6.2830 5 
= #8 
5 6 
Const Solution 
n Cr 
- ,219 086 
2,439 .074 
-.658 $118 
2,987 ,166 
-1,516 ¿20€ 
-1,974 2990 
-2,631 482 
-3,288 610 
-,110 022 
-,219 .045 
- ,529 „ОТО 
2,494 110 
2,658 . 151 
-.987 250 
-1.316 +960 
-1,644 „480 
= „О55. 012 
-,110 2026 
- „164 ‚040 
- ,247 .064 
-.929 ,090 
2,494 2195 
- .658 „287 
2,822 2590 
2,087 .009 
-.073 „О17 
=„110 .028 
- .164 „044 
-„219 .060 
-.329 LO 
-.238 170 
- ,548 „240 





= 


10 


TABLE [TÍ Cont 'd.) 


Š 


Coeff of 


C. e 


P 


-2 .431 
-2 .358 
-2 „285 
-2.175 
-2 065 
-1,845 
-1.626 
-1,407 


-2 „404 
-2.950 
-2,255 
-2,144 
-2 ,052 
-1.808 
-1,585 
-1.562 


4 


Coeff of 


it 


q 


4.118 
$4985 
3.852 
5.656 
5 .454 
5.056 
2.659 
СИОТ 


4 ¿062 
5.924 
5.786 
5.079 
5.572 
2.957 
2.543 
22123 


5 


Const 


—„027 
- ¿059 
- ¿082 
-.125 
-.] 64 
-.247 
= 069 
- 411 


=„О22 
- „044 
- ¿066 
= ¿099 
-= 0.132 
-,197 
- 263 
=. 529 
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Solution 


¿007 
2014 
Ос 
0035 
.049 
.085 
.154 
. 200 


.005 
ol 
.018 
-029 
040 
-069 
„110 
„160 
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TABLE V 


Wave Rise "Correction" Factors for Wetted Length Determinations 


from Heave Readings 


Angle Photo Velocity Speed 


of 
Trim 
> 


(дер .) 


18° 


21° 


30 


No, 


V 
(ft/sec. ) 


27.72 
27.72 
27.78 
27.78 
21.12 
2 
27.78 
27.78 


27.72 
27.72 
21.12 
18 
27.66 
27.12 
27.712 
27.12 


27.72 
27.84 
27.72 
27 (2 
27.72 
27.72 
27.90 


11.96 
11.96 
11.99 
11,99 
11.96 
11.96 
11.99 
11.99 


11.96 
11.96 
11.96 
11.96 
11.9) 
11.96 
11.96 
11.96 


11.96 
12.02 
11.96 
11.98 
11736 
11.96 
12. ОВ 


Wetted 
Coeff. Length/Beam 


Ratio 


from Photo 
А. 


2.30 
2. 30 
SO í 
5.05 
6.55 
6.13 
1822 
Telih 


2.68 
2.63 
5.87 
5.82 
7.25 
7.12 
3.80 
3.80 


2.52 
2.12 
5.72 
5.67 
1.03 
So 
1.00 
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TABLE VI 


Reduction of Experimental Data, Correlating Cr and À 


$ = .97 slugs 
24.3 

Pict V 
519 27.72 
550 27.72 
551 27.78 
552 27.78 
53 27.72 
554 Ale? 
555 27.78 
556 27.78 
557 27.712 
558 27.72 
559 27472 
560 27.12 
561 27.66 
562 дета 
563 27.72 
56h 2 
565 27.72 
566 27.84 
567 DUO 
568 27.72 
569 27.72 
570 21212 


27.90 


b 


2 





745.35 
715.35 


715.35 
751.82 
715. 35 
745.35 
715.35 
715.35 
755.06 


02777 ғ. 


À 


o 


9 


мм до O 


o 


— 3 O. ONN IV, № 
e 
Б о рл оом м 


o 
со CO ON ON 
мо 


=з \л ул. гә N 
o 

N 

VA 


T2 
3,80 
3,80 


2.52 


1.03 
li. OO 


q w 


ea], 


17.60 


55 


17.6055 


105. О 
10h. 95 
135.58 
133112 
150 ° 09 
151.66 


55.47 
Sho lhi 
121.19 
120.115 
1119. 39 
117. 36 
78.63 
78.63 


52.16 
50.52 
118.36 
117.39 
145.49 
131.85 
83.88 


1..2 В Tt. 
56! 2 


10.61 
10.61 
18.05 
18.05 
21.8) 
21.81, 
2l. 71 
271 


1l. 81 
1}. 81 
27.54 
27.54 
32.76 
32.76 
19.93 
19.93 


18.93 
18.93 
36.69 
36.69 
11. 10 
0.28 
27.13 


.223 
.223 
ci 
„172? 
‚161 
. 16h 


.160 


.267 
.272 
.227 
.229 
„219 
„222 
.253 
253 


.363 
. 375 
„310 
.313 
» 303 
o 305 
2337 


бб, 





TABLE VII 


Tabulation of Computed and Experimental Results for 
the Gravity Independent Lift Coefficient 


1 2 Š 4 5 6 d 8 


Definition of 
ñ Based on X, 


Definition of 
ñ Based on Ç 


^L Comp “L Comp 
BR Z C C C C C C 
V L Exp L Comp 1 Exp 1 Сотр L Exp 
ES 2° 21.44 „051 Oar аиб „ОЛ? 2052 
„40 18.51 031 „060. 1,935 .065 2.10 
45 24,61 ,031 ‚058 1,87 бы 205 
2 20.53 ‚035 „049 1.485 „052 1.58 
„68 20.74 „029 ‚046 1.585 ¿049 1.69 
„70 24 „92 „029 (045 1 ,048 1.65 
.87 14,55 ‚023 ‚041 1,78 „041 1.78 
„87 22.88 ,028 ‚oa T6 „041 1.46 
217 20.04 „027 035 1.30 ‚033 122 
1.57 25.01 „022 029 1.32 ‚027 1,23 
2.54 13.37 „019 „022 1,16 020 1,05 
3 00 13.66 „015 ¿020 1,33 ,018 1920 
5.00 1з 7° „015 020 11% „018 1820 
5 .00 16.99 „015 „020 1.33 018 1.20 
2.00 21,87 „015 „020 1.33 „018 1.20 
4.79 12.20 „012 ‚015 1,25 „012 1,00 
6.50 25.01 .009 Olle 1.33 ‚008 .89 
6.50 18.64 „010 Ole J OO .008 ,80 
7.79 13.39 ,009 „010 1.11 .005 „555 
8.12 23 „39 „009 „010 1.11 „005 „555 
8„27 12 „78 ‚009 „010 J 11 ,005 „555 
Ес 14,40 ,082 о ea ‚150 1.83 
,85 19,98 ,091 „125 1287 as пи: 
„38 16.99 117 „125 1.05 Glas 1.14 
„58 24 „86 „091 „125 1.35 133 1.46 
„50 13.57 „092 „120 1. 120 1.3 
„54 13 „68 „084 „108 ies ,115 137 
62 21956 ‚073 100 er 105 za 
„85 24 „92 073 ,085 1,16 „088 1.20 
1.45 21.87 .055 063 1,14 ‚057 1,04 
1,45 12.81 .054 ‚063 TTE, ,057 1.06 
1.50 le e „054 ‚062 IL us os 104 
1.52 16.23 „053 — NE ‚055 1,04 
1.64 ler „048 ‚059 1,23 ‚053 128 
№ 66 12,86 ,049 059 1,20 „052 1,06 
2.20 24 ,80 ‚041 бз aa „042 1,02 
3 „66 10.97 029 039 1,34 ,029 1,00 
4.07 p ,029 Оо i 1026 „895 
4,54 13,48 „026 „054 1.31 „024 „925 
з 20,95 „027 „054 1,26 „025 „854 
4.66 18.13 „025 „055 1„32 ‚023 „920 
4 ,82 151 „026 055 J 7 „022 „846 
5,14 22.97 ‚027 „052 1.18 ‚021 „778 
5,54 17,57 ‚022 ios] E ‚020 ‚910 
(Т 16 12 „00 021 029 1.38 „020 „955 
8,54 15,46 ,019 1025 1.37 О „650 
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O 


195,89 
24,77 
16,68 
24 ¿92 
12.69 
16.32 
212.90 
ей 
25.10 
12.96 
24,46 
24 .10 
17.60 
20.86 
1797 
17.60 
10.85 
24 „16 
14.55 
19.86 
ar 
10,72 
24,46 
20.74 
17.08 
12,41 
15701 
17.51 
19290 
1627 
14,13 


15.70 
11,89 
18,28 
2250 
14.10 
15.40 
15.78 
11„98 
10982 
12.50 


24 ,98 
20.04 
16.96 
24 .95 
2159 


TABLE VII(Cont 'd, ) 


4 5 


6 


Definition of 
m Based on T 


L Exp Sp Comp 
. 147 2207 
259 2205 
‚ 15$ . 194 
‚165 JE 
143 . 158 
„145 . 15838 
ДОО o 150 
o122 149 
ПЁ] °148 
2096 En 
2086 2109 
088 o 104 
JOY , 104 
084 , 100 
‚076 „095 
‚060 2075 
‚059 ¿Oe 
057 ‚070 
‚055 069 
‚055 ‚069 
052 068 
-029 ‚065 
‚045 ‚056 
044 ‚055 
044 095 
044 0055 
‚045 ‚054 
‚059 ‚055 
0037 051 
‚056 ‚051 
0056 ‚050 
о 19% 
leg 2197 
ЖТ „146 
001 , lel 
090 , 120 
084 „114 
„О7О o 104 
о o 104 
‚066 ‚096 
‚066 ‚096 
‚256 ‚455 
2518 2445 
2971 ‚455 
‚285 0435 
905 ‚425 


C 
L Comp 


E Exp 


1,41 
1,47 
121 
111 
1.10 
1.03 
1922 
1222 
1.15 
1.12 
ieee 
1218 
1.20 
115 
1.25 
1.22 
1522 
1.25 
1525 
1225 
Ix 
1.33 
1,24 
1.25 
125 
1.25 
1,26 
1,36 
1.58 
1,42 
1:99 


1.15 
1.18 
1.25 
1.55 
1.55 
1.56 
1.49 
1.47 
1.45 
1,45 


1,78 
1.40 
І 7 
1,55 
1259 


7 


8 


Definition of 
m Based on e, 


.208 
205 
199 
2184 
1956 
2156 
„146 
„145 
,144 
0098 
094 
0093 
‚095 
¿0390 
.085 
053 
‚052 
‚050 
049 
048 
.047 
‚045 
‚ОТ 
0030 
029 
‚029 
028 
‚025 
O21 
0020 
„О2О 


,180 
» 180 
147 
OBI 
.080 
2069 
2048 
048 
057 
‚056 


,4 60 
445 
2455 
0455 
2425 


4 Comp 
C C 
L Comp `L Exp 


1.42 
1.48 
127 
1215 
1.09 
1,08 
1,19 
12219 
1.10 
1,02 
1592 
1.06 
107 
1.07 
Im 
084 
„880 
808 
09 
.87 
o YO 
.88 
299 
.67 
.66 
.66 
o. 69 
„64 
ый 
0955 
0955 


1.05 


1 17 
1,53 
1,39 
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TABLE VII (Cont 'd.) 


5 4 9 6 T 8 


Definition of 
m Based on Ç 


Definition of 
m Based on A 


12 


18 


: a en Comp Cd, Comp 
V L Exp d: Comp L Exp %, Сотр Е Ехр 
29-16 .302 2400 1.53 2410 1.56 
129375 393 2400 1,02 410 1.04 
oma ‚o24 sooo 1.19 930 1120 
21259 „294 2580 1. 0085 1.31 
28.25 29] 20 10 1927 375 1429 
24 „98 . 198 Sort 1285 sone 1.88 
13.79 „280 2945 1.48 2345 Jn 5 
Т „дё . 918 2345 1.08 345 1.08 
21.14 216 29595 121 385 1.21 
24,40 .254 928 1.2 .350 1559 
12228 268 919 Т 19110 1418 
J2 22 „264 . 918 1.40 ¿One 1218 
17,54 at 2a als 2175 ON 
24,40 2258 „288 21 268 Те 
20.07 2240 2208 1,20 ‚268 1211 
10 ,80 2695 a205 121 „266 1.15 
24 ,49 ,204 2248 1:22 2214 1.05 
20.89 2209 2247 1.15 als J.G: 
14,95 ¿206 2245 1,18 „210 1„02 
12.38 „204 ‚242 1.29 .208 1202 
12,50 .195 ,240 1.25 ‚205 1.05 
24 „28 o 164 2210 1229 . 164 1.00 
асет ler „204 1827 155 294 
14,94 2157 ¿206 1929 150 „96 
M „154 -201 1.51 „149 ЗЭ 
16.78 „ 145 O 187 JS „128 87 
12922 „140 „185 152 sled 299 
ЕО 156 «EOM 1 55 116 .85 
11.74 188 „180 1551 15 .83 
20.74 .136 „180 1,92 114 284 
16,23 o 134 1279 1.54 1.15 „84 
116 505 adole SH 1,35 El 05 
T6978 125 oom 1256 ‚088 2 
1052 022 Gy 1097 .087 Tl 
JET „КО 1 ,166 1257 .086 SUM 
14,55 ,120 1658 И ‚085 211 
202 Ө ¿109 1,40 .085 e14 
25.10 2509 .585 1545 . 980 1,14 
1263 441 9 60 ] e 7 . 960 1227 
zu S 2450 2560 1,24 » 060 1.24 
25.50 ,4 78 ,960 1,18 . 960 1 18 
10.15 „622 550 ‚885 0900 88 
21,26 2490 oS 1.09 EU 1.07 
15,49 2342 29595 1,56 a923 1.55 
10.13 444 25219 I Eo >15 1.16 
1728 2461 92/9 1,14 solo Jall 
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21.01 
24,43 
24,40 
17.68 
10,74 
1083 
19,82 
15.91 
24,61 
20.83 
12,44 
25.18 
10.98 
14,91 
2.01 
16.51 
20.98 
16590 
15,65 
12,09 
16.18 
14.61 
14.70 
18,64 
16220 
11,90 
11,30 
11290 
113799 
11,99 
11296 
965 
11.99 
11298 


gli lil 
17.85 
13,02 
10 „10 
17,45 
20 ¿04 
24 40 
24 „34 
10,80 
24.61 
20.86 
12,47 
SUI 


TABLE VII (Cont 'd.) 


4 5 


6 


Definition of 
m Based on Ç 


0 

L Comp 

E Exp on Comp i Exp 
„448 ‚529 1.18 
0433 1529 1822 

od 005 0867 
„446 „492 1,10 
„442 „492 1. №1 
‚432 488 1.1% 
„419 „482 1.15 
424 ‚456 1208 
«400 „456 1„14 
„409 450 LO 
„581 441 12276 
s961 408 1.15 
946 406 J1 
03843 „405 МИ 
some 695 1225 
.296 .968 1526 
„2817 566 1327 
.276 2562 1.91 
„290 „559 1.24 
.276 SOS 1.25 
.266 ‚943 1249 
.265 3495 1,29 
„249 009 1296 
„249 вой J 55 
„247 odo 1.56 
1255 od UI 1,53 
220 .328 1,47 
¿225 „528 1.47 
171 „284' 1„66 
„Lie „284 1265 
„164 21.9 1958 
.161 .274 1.70 
2:165 sa te 1:65 
,160 оТ 1289 
.565 „680 1.20 
.579 . 680 1. №5 
923 „бё 1729 
2501] ¿671 1,34 
.568 ,645 1,14 
¿2000 „645 115 
.558 .645 1516 
„473 „боб 1.54 
480 695 192 
„527 „618 117 
„544 .617 1.15 
. 995 610 1.21 
„492 . 605 1223 


7 


8 


Definition of 
m Based on o 


‚515 
.915 
487 
.470 
470 
2465 
„460 
„426 
„425 
„420 
„416 
2068 
964 
‚961 
29590 
оО 
907 
3299 
.294 
DoD 
„261 
2600 
2254 
2250 
2250 
2250 
2280 
2230 
12 
. 128 
2104 
ТОЕ 
094 
s091 


.686 
. 686 
„675 
2675 
o 650 
2650 
. 650 
«боб 
2686 
-615 
,614 
,608 
„604 


Gi Comp 


C C 
L Comp L Exp 


1,15 
1229 

2905 
1.05 
1206 
1207 
1210 
1,01 
1,06 
1.03 
1709 
1,02 
1.05 
1.05 
lowe 
1,05 
1.07 
1.08 
1201 
1.05 

„98 

auc 
12092 
1.00 
101 

-99 
1.05 
1.08 

о 14 

74 

. 65 

„бә 

sot 

ool 


1.21 
1:19 
1:29 
1.55 
1.15 
1.15 
Ir 
1.55 
1,33 
1.17 
1215 
121 
1.23 


ox. 





TABLE VII(Cont'd,) 


ó 4 5 6 T 8 


Definition of 
th Based on © 


Definition of 
m Based on A 


24 


30 


C б 
à m Comp L Comp 

Cv 1 Exp L Comp L Exp a Comp Cr Exp 
24,95 ,510 ,594 1.16 .585 1.15 
2077 „448 ,594 1,32 ,585 1.30 
12,41 ,453 ,594 1591 ,585 ско! 
18,39 ,484 ,578 1.20 . 560 1.16 
14,88 .467 „570 1.22 ,550 1.08 
20,92 ,459 „570 1.24 „550 1„20 
23005 ‚458 „570 1„25 „550 1„20 
11,65 „454 „559 1„24 . 500 1,15 
13,42 ,423 .539 1.27 ,498 me: 
20,62 „419 „5856 1„28 „495 1018 
16,16 ‚408 GEE 1.31 „487 1,19 
16.59 ‚437 ‚517 1,18 457 1,05 
18,54 ,416 ,514 1,23 ‚450 1.08 
11,90 „419 „514 1„25 „450 1„07 

„599 „510 1.28 „445% nes 
11,96 „272. „452 1.66 .308 1.13 
11,96 „267 „451 1,59 ‚306 1535 
11,96 „255 „450 1.70 ,245 „97 
Tess „229 AN 1,80 ,175 Ub 
11,96 .227 sae 1.81 .174 ST 
11,96 ‚222 „406 1S ,148 ,67 
11,94 „219 „405 1„85 „146 „67 
2172 „759 „826 1212 ‚835 1 15 
16.41 79 ,826 1.15 .835 1.16 
12.78 Tid „826 1.16 .835 TNNT 
24,98 „740 „822. 1211 .830 1.12: 
20,86 „796 „814 1,02 .820 1.05 
17,63 Brora ,814 1,05 ,820 1,06 
10,13 „692 ,807 1917 ‚816 1,18 
17,45 ,674 „783 1,16 „785 1.16 
10,89 ,665 „78% 1.18 ,785 1.18 
20 ,10 .578 „770 ПЕ „765 ee 
12,41 ,622 „749 1 20 „749 1.20 
14,91 .623 „749 1020 „749 1,20 
14 ,98 .598 „718 1„20 SION 1.19 
13,24 ,566 ,686 121 „668 1,18 
11,99 „514 „658 1.28 ,622 1427 
14 „49 „518 „658 DEDI ,622 1.20 
12 ,02 0375 .583 1.55 ‚443 1.18 
11,96 5065 „581 1„60 „455 1.20 
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Values of the Lift Coefficient for Low Aspect Ratio Airfoils 
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Schematic for Spray Thickness, 
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Schematic for Heave 
Calibration Dimensions 


Fig. £_ 





Undisturbed 
Actual Water Surface 


eeu р Р Н эз 





Lift of a Two- Dimensional Flat 


Plate as a Function of Froude Number 
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Figures for APPENDIX II 
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